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Abstract. This paper presents a fractal analysis of radiographic patterns obtained from specimens 
with three types of inserted welding defects: lack of fusion, lack of penetration, and porosity. The 
study focused on patterns of carbon steel beads from radiographs of the International Institute of 
Welding (IIW). The radiographs were scanned using a greyscale with 256 levels, and the fractal 
features of the surfaces constructed from the radiographic images were characterized by means of 
Hurst, detrended-fluctuation, and minimal-cover analyses. A Karhunen-Loeve transformation was 
then used to classify the curves obtained from the fractal analyses of the various images, and a study 
of the classification errors was performed. The obtained results indicate that fractal analyses can be 
an effective additional tool for pattern recognition of weld defects in radiographic tests. 

Keywords: Welding defects, pattern recognition, fractal analyses 
PACS: 42.30.Sy,81.20.Vj,87.59.Bh,87.59.Hp 

1. INTRODUCTION 

In a recent paper, Silva et al. [1] discussed the characterization of failure mechanisms 
that occur in fiberglass-reinforced polymeric matrix composites when subjected to ten- 
sile and flexural loads. The characterization was based on the analysis of acoustic- 
emission signals emitted by the composite during the process of failure. By looking 
at some fractal properties of the acoustic emission signals, namely, the fractal indices 
related to the Hurst analysis detrended-fluctuation analysis [3], minimal cover anal- 
ysis [0], and the box-counting dimension analysis [5], they were able to distinguish the 
different failure modes. 

The study presented in this paper aims to characterize, through fractal analyses, the 
welding defects present in radiographic images. We focused on patterns of carbon steel 
beads from radiographs of the International Institute of Welding (IIW). The images were 
scanned in 8-bit resolution (256 levels of grey), and then processed by using the software 
Image Pro Plus 4.0. In order to improve contrast, a median type filter was used to smooth 
unpredicted noise. From the scanned images, we built surfaces by associating the grey 
level at each pixel with a height variable. The patterns in each surface were then studied 
by fractal analyses. 

In order to establish the parameters to be calculated, we first present a brief review 
of the numerical analysis used in the treatment of the data. Afterwards, we present and 
discuss the results obtained. 
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FIGURE 1. The bottom shows a radiographic image of a defect (lack of fusion), and above it is the 
corresponding surface, obtained by converting the greyscale to height variables. 



2. FRACTAL ANALYSES 

The numerical treatment of the images was performed on data generated from the 8-bit 
scanning of the radiographs, corresponding to 256 levels of grey, which are translated 
into a height variable Zij. Here, represents the coordinates of a pixel, with i = 1, 
2, L x and j = 1, 2, L y for an image containing L x x L y pixels. This process is 
illustrated in Fig. Q] 

In the fractal analyses we considered extended two-dimensional versions of the Hurst 
(or R/S) analysis [2], detrended-fTuctuation analysis (DFA) and minimal-cover anal- 
ysis [0]. In order to make the paper self-contained, we will present the details of the 
numerical techniques used in the analysis of the surfaces associated with the images, but 
first we make a few remarks that are generally valid. 

All techniques start by covering the image with a grid composed of square cells 
containing txt pixels, making sure that the grid is centered in both directions. This 
guarantees that, if L x or L y are not multiples of t, only pixels in the periphery of the 
image are left outside the grid. Each technique then involves the calculation of the 
average of some quantity Q{t) over all cells, for different values of T. In a surface with 
genuine fractal features, Q(x) should scale as a power of x for T> 1, 

Q(x)~x n - 

Fractals of different nature should give rise to different exponents rj, providing a sig- 
nature of the fractal. In our case, due to the finite amount of pixels and the limited 
resolution of the heights ztj, this power-law behavior is hard to observe. However, as 



discussed in the final section, we can still use the T-dependence of the functions Q(l) to 
characterize the different defects. 



2.1. Hurst analysis 

The rescaled-range (R/S) analysis was introduced by Hurst |0] as a tool for evaluating 
the persistency or antipersistency of a time series. The method works by dividing the 
series into a series of intervals, and calculating the average ratio of the range (the 
difference between the maximum and minimum values of the series) to the standard 
deviation inside each interval. The size of each interval is then varied. 

We define a two-dimensional version of the R/S analysis in the following way. Given 
a T x t cell, whose lower left corner is located at pixel (io,jo), we calculate (z) T , the 
average of Zij inside the cell, 

2 y 'o+t-Uo+t-I 
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where the primed summation runs over all pixels (i, j) inside the cell. We then define an 
accumulated deviation from the mean as 

Zij= E E (*«-<*>*). (2) 
k=k l=jo 

from which we extract a range, 

R(t) = max Z ;J - — min Z ; y (3) 

;'o<i</o+T-l ('o</^;'o+T— 1 

and the corresponding standard deviation, 



W = ^E'4 w 

V (i-J) 

Finally, we obtain the rescaled range R(t)/S(t), and take its average over all cells. 

In a surface with true fractal features, the rescaled range should satisfy the scaling 
form 

where H is the Hurst exponent. 

A typical curve obtained from the R/S analysis of the surfaces is shown in Fig. [2a). 



2.2. Detrended-fluctuation analysis 



The detrended-fluctuation analysis (DFA) aims to improve the evaluation of cor- 
relations in a time series by eliminating trends in the data. 
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FIGURE 2. Typical curves obtained from the fractal analyses. The quantity L is defined as L — yjL x L y . 



Our two-dimensional extension of the method consists initially in obtaining a new 
integrated two-dimensional data set Zij, 

fy = ££(*«-<*»> (6) 

where the average (z) is taken over all pixels, 

After building the grid (with cells of side t), the integrated data inside a given cell is 
fitted by a plane. Then, a detrended variation function A ;; - is obtained by subtracting 
from the integrated data the local trend as given by the fit. Explicitly, we define 

A ij=Zij-hij, (8) 

where htj is the height associated with pixel according to the fit. Finally, we 
calculate the root mean-square fluctuation F(t) inside a cell as 



V iu) 

and average over all cells. For a true fractal surface, F(t) should behave as 

F(T)~T a , (10) 

where a is the scaling exponent. 

A typical curve obtained from the detrented-fluctuation analysis of the surfaces is 
shown in Fig.|2tb). 



2.3. Minimal-cover analysis 



This method has been recently introduced [4], and, in its original version, it relates the 
minimal area necessary to cover a given plane curve, at a specified scale, to a power-law 
behavior. 

In our two-dimensional extension, the method relates the minimal volume necessary 
to cover a given surface, at a specified scale. After building the grid, we can associate 
with each txt cell, labeled by a variable k, a prism of height A^, defined as the 
difference between the maximum and minimum values of zu inside cell k, 

A k = max zij- min Zij. (11) 

The minimal volume is then given by 

V(t)=J>t 2 , (12) 

k 

where the summation runs over all cells. 

Ideally, in the scaling region (t>1),V(t) should behave as 

V(t)~t 3 ^, (13) 

where is the minimal cover dimension, which is equal to 2 when the surface presents 
no fractality. 

A typical curve obtained from the minimal-cover analysis of the surfaces is shown in 
Fig. He). 



3. RESULTS AND DISCUSSION 

In order to classify the images, we used a supervised variation of the Karhunen-Loeve 
(KL) transformation @], and applied it to the set of curves produced by the fractal 
analyses described in the previous Section. In this sense, the fractal analyses can be seen 
as a sophisticated preprocessing tool. For each image, we collected the corresponding 
curves from the three fractal analyses, forming a single vector with M components 
(AT = 57). 

We proceeded by first randomly dividing the vectors into a training set (with N = 120 
vectors) and a test set (with 31 vectors), performing all the relevant operations (as 
described below), and calculating the confusion tables. Finally, we took averages over 
300 different choices of training and test sets. 

Let x,- be the vector corresponding to the ith image. The KL transformation consists of 
first projecting the training vectors along the eigenvectors of the within-class covariance 
matrix Sty, defined by 

I N C N k 

Sw = m ^ Y,yik( x i- m k)( x i- m k) T : (14) 



TABLE 1. Average confusion matrix for the training 
vectors. The possible classes are lack of fusion (LF), lack 
of penetration (LP), porosity (PO) and no defects (ND). 
The figures in parenthesis indicate the standard deviations, 
calculated over 300 sets. The value in row i, column j 
indicates the percentage of vectors belonging to class i 
which were associated with class j. 





LF 


LP 


PO 


ND 


LF 


87.4 (4.8) 


7.5 (3.6) 


4.9 (4.3) 


0.2 (0.9) 


LP 


1.5(2.1) 


93.8 (3.4) 


2.1 (2.6) 


2.5 (2.4) 


PO 


0.2(1.0) 


1.4 (2.7) 


97.0 (3.8) 


1.4 (2.6) 


ND 


1.3(1.2) 


3.7(1.8) 


1.6(1.1) 


93.4 (2.3) 



where Nc = 4 is the number of different classes, is the number of vectors in class 
k, iti£ is the average vector of class k, and T denotes the transpose of a matrix (in this 
case, of a column vector). The element is equal to one if x,- belongs to class k, and 
zero otherwise. We also rescale the resulting vectors by a diagonal matrix built from the 
eigenvalues A; of Siy . In matrix notation, this operation can be written as 

X' = A-2U r X, (15) 

where X is the matrix whose columns are the training vectors x,-, A = diag(Ai, A2, ...), 
and U is the matrix whose columns are the eigenvectors of S\v- This choice of coordi- 
nates makes sure that the transformed within-class covariance matrix corresponds to the 
unit matrix. Finally, in order to compress the class information, we project the resulting 
vectors onto the eigenvectors of the between-class covariance matrix Sg, 

Sb= E TT-f^-mXm^-m) 7 , (16) 
k=\ ;v 

where m is the overall average vector. The full transformation can be written as 

X" = V r A~3U r X, (17) 

where V is the matrix whose columns are the eigenvectors of Sg (calculated from X'). 

With 4 possible classes, the fully-transformed vectors have 4 — 1=3 relevant compo- 
nents [6]. A vector x,- is associated with the class whose average vector lies closer to x, 
within the transformed three-dimensional space. 

After using the training vectors to obtain the full transformation, the classification 
can be checked with those same vectors. The resulting average confusion matrix is dis- 
played in Table[TJ When employed to classify the testing vectors, the KL transformation 
produces the confusion matrix shown in Table [2] From the tables we can see that the 
training vectors are associated with the correct class in about 90% of the cases, whereas 
the percentage of correct classification for the testing vectors is about 50%. This last 
result surely lies above the expected rate produced by random association (which would 
be 25% in this case). 



TABLE 2. The same as in TableQ] for the testing vectors. 





LF 


LP 


PO 


ND 


LF 


48.1 (23.6) 


19.8 (20.7) 


22.8(19.9) 


9.2(14.7) 


LP 


10.7(11.9) 


55.1 (18.8) 


14.3 (14.0) 


19.9(16.1) 


PO 


12.2(16.3) 


14.7(18.6) 


57.6 (25.0) 


15.5(19.1) 


ND 


9.4 (8.7) 


19.2(11.2) 


13.7(10.1) 


57.6(14.3) 



In summary, we have shown that fractal analyses are a promising tool for classifying 
welding defects in radiographic images. We believe that its efficiency can be consider- 
ably improved by using a larger sample of images, as well as by working with 16-bit 
scans, which would greatly enhance the grey-level resolution. We are currently working 
on this direction. 
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